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We show that the hnearly polarized gluon distributions appear in the color dipole model as we 
derive the full cross sections of the DIS dijet production and the Drell-Yan dijet (7* jet correlation) 
process. Together with the normal Weizsacker- Williams gluon distribution, the linearly polarized 
one will contribute to the DIS dijet production cross section as the coefficient of the cos (2A(/!>) term 
in the correlation limit. We also derive the exact results for the cross section of the Drell-Yan 
dijet process, and find that the linearly polarized dipole gluon distribution which is identical to the 
normal dipole gluon distribution involves in the cross section. The results obtained in this paper 
agree with the previous transverse momentum dependent factorization study. We further derive the 
small-s evolution of these linearly polarized gluon distributions and find that they rise as x gets 
Ch ■ small at high energy. 
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I. INTRODUCTION 

Q-1 In small- a; physics, two different unintegrated gluon distributions [l|, 0] (also known as transverse momentum 
pi I dependent gluon distributions), namely the Weizsacker- Williams gluon distribution xG^^^ [SHi and the dipole gluon 
. distribution xG*^^-* , have been widely used in the literature. The Weizsacker- Williams gluon distribution can be 
• interpreted as the number density of gluons inside dense hadrons in light-cone gauge. The dipole gluon distribution, 
] despite of lacking the probabilistic interpretation, has been thoroughly studied since it appears in many physical 
. processes 0, Q and it is defined via the Fourier transform of the simple color dipole amplitude. This dipole gluon 
^ ' distribution can be probed directly in photon-jet correlations and Drell-Yan dijet measurement in pA collisions, 
on . Recent studies 0, Q on the Weizsacker- Williams gluon distribution indicate that it can be directly measured in DIS 
0^ ■ dijet production and its operator definition is related to color quadrupoles instead of normal color dipoles. Other 
' more complicated dijet processes in pA collisions (e.g., qg or gg dijets) involve both of these gluon distributions 
^\ • through convolution in transverse momentum space. The complete calculations were performed in Ref. Q in 
On [ both the Transverse Momentum Dependent (TMD) factorization formalism and the color dipole model. The results 
■ demonstrate the complete agreement between these two formalisms in the kinematical region where they are both 
valid. 

Linearly polarized gluon distributions, denoted as h^]^ {x, q±'j^, where x and q± are the active gluon's longitudinal 
momentum fraction and its transverse momentum, respectively, were first introduced in Ref. [9]. This new gluon 
distribution effectively measures an averaged quantum interference between a scattering amplitude with an active 
gluon polarized along the a:;(or ?/)-axis and a complex conjugate amplitude with an active gluon polarized along 
[ the y{oT x)-axis inside an unpolarized hadron. Because of the unique transverse spin correlation between the two 
gluon fields of the distribution, the linearly polarized gluon distribution can contribute to a physical observable with 
cos(2A0)-type azimuthal angular dependence, or the azimuthally symmetric observables if they come in pairs. As 
proposed in Ref. [loj . the linearly polarized gluon distributions can be directly probed in dijet and heavy quark 
pair production processes in electron-hadron collisions. As expected, this distribution also contributes to the cross 
section in photon pair productions [ll|, IHl and the Standard Model Higgs boson production p^l - fTsj in hadron- hadron 
collisions. Since the integrated parton distributions for incoming protons were used in the calculations of pA collisions 
in Ref. the linearly polarized gluon distribution does not enter the cross section except for the Drell-Yan dijet 

processes as we show in the later discussion. 

In Ref. [1^, the linearly polarized partner of both Weizsacker- Williams and dipole gluon distributions inside an 
unpolarized nucleus target is studied in Color Glass Condensate (CGC) formalism. The corresponding cross sections 
of deep inelastic scattering (DIS) dijet production and the Drell-Yan processes in pA collisions are computed in terms 
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Normally it is denoted as ^{x,qj_). Here throughout the paper, in order to avoid confusion on the notation, we use {x, qj_) with 
i = 1, 2 to represent the hnearly polarized gluon distributions. 
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of the TMD formahsm. In both processes, the hnearly polarized gluon distributions appear as the coefRcients of the 
cos(2A0) term in the cross section, where they were found to be consistent with the small- a; formalism as well [Tcj . 

Inspired by Ref. [ll], we perform the detailed calculation in the color dipole model for the DIS dijet production 
and the Drell-Yan dijet processes in pA collisions and we find identical results as those with the TMD formalism 
for the cross sections in the correlation limit, which is defined as a limit when the final state dijets are almost back- 
to-back. For the DIS dijet production, the complete analysis of the quadrupole amplitude shows that the linearly 
polarized gluon distribution of the Weizsacker- Williams type comes from the off-diagonal expansion of the quadrupole 
amplitude. Using a hybrid factorization, we obtain the exact results for the cross section of the Drell-Yan processes 
in pA collisions. In the correlation limit, this exact result reduces to the TMD cross section obtained in Ref. [l6| . 

Another objective of this paper is to study the small-x evolution of the linearly polarized gluon distributions. 
The small- a; evolution of the dipole type linearly polarized gluon distribution is essentially the evolution of the dipole 
amplitude, which is governed by the Balitsky-Kovchegov equation Derived from the evolution of quadrupoles, 

the evolution of the linearly polarized Weizsacker- Williams gluon distribution is quite complicated. Nevertheless, in 
the dilute regime, we find that both linearly polarized gluon distributions receive the exponential enhancement in 
terms of rapidity at high energy as the normal unpolarized gluon distributions do due to the small- a; evolution. 

The rest of the paper is organized as follows. In Sec. II, we calculate the cross sections of the DIS dijet production 
and the Drell-Yan processes in pA collisions and demonstrate that the linearly polarized gluon distributions naturally 
arise in the dipole model. We discuss the small- a; evolution equations of the linearly polarized gluon distributions in 
Sec. III. The summary and further discussions are given in Sec. IV. 



II. THE LINEARLY POLARIZED GLUON DISTRIBUTION IN DIPOLE MODEL 

In this section, following Ref. Q, we show that the cross section of the DIS dijet production and the Drell-Yan 
dijet process in the color dipole model, namely the CGC approach, involves the linearly polarized gluon distribution 
as well. The reason why this does not appear in the original work in 's'l is that there the azimuthal orientation of the 
outgoing partons was averaged over. 



A. DIS dijet production 



After averaging over the photon's polarization and summing over the quark and antiquark helicities and colors, the 
cross section of the DIS dijet production in the color dipole model can be cast into 



d^xi d^x'i d^X2 d^x'2 



(27r)2 (27r)2 {2tt)^ (27r)2 



XafI 



1 + , 2:2 ; ^2 , a;'i ) - (xi , X2 ) - S^^ (x'^ , x[ ) 



(1) 



where the two- and four-point functions, which are characterized by the Wilson lines, take care of the multiple 
scatterings between the qq-pair and the target. They are defined as 



S^^J{xi,X2) - ^ (TrC/(xi)C/t(x2))^^ , 
5W(xi,a:2;x^,x;) = ^^{TtU{x,)U\x[)U{x',)U\x2))^^ 
with U{x) ^Ve^j) l^igs J da;+ T'^ (a;+ , a;) | . 



(2) 
(3) 
(4) 



The notation {. ■.)x is used for the CGC average of the color charges over the nuclear wave function where Xg is 
the smallest fraction of longitudinal momentum probed, and is determined by the kinematics. The splitting wave 
function of the virtual photon with longitudinal momentum p+ and virtuality in transverse coordinate space takes 
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the form, 



«e/^i(e/kl)-|7r['^"+'^'3+(l - z)^ (5„_(5^_2], A = 1, 
ze/i^i(e/|r|)^[,5o_,50_(l ~ z) + <5„+<5^+0], A = 2, 



V'^^(p+,z,r) = 2^^/— z(l - z)Qifo(e/|r|)(Sa;3. 



(5) 



(6) 



where z is the momentum fraction of the photon carried by the quark, A is the photon polarization, a and /3 are the 



quark and antiquark hehcities, r the transverse separation of the pair, e? — z(l — z)Q^, and the quarks are assumed 



to be massless. 

In order to take the correlation limit, we introduce the transverse coordinate variables: u — xi — X2 and v = 
zxi + (1 — z)x2, and similarly for the primed coordinates, with respective conjugate momenta P± = (1 — z)ki± — zfc2jJEI 
and q±. The correlation limit {P± ~ ki± ~ /c2_l 3> q±) is therefore enforced by assuming u and u' are small as 
compared to v and v' and then expanding the integrand with respect to these two variables before performing the 
Fourier transform. Following the derivation in Ref. [sj, one can find that the lowest order expansion in u and u' of 
the last line of Eq. ^ gives 



With the help of the identities 



d'^u cPu' 



1 



(27r)2 (27r)2 



Uiu'^Ko{efu)Ko{€fu') = 



(2^)2 
1 



2\2 



_iPl+ej) 



(Pi +4)4 



(2^)2 (27r)2 'J ^ ' ^ ' (2^)2 (P2 +g2)4 ' 

one can integrate over u and u' and obtain the complete differential cross section in the correlation limit. 



(7) 

(8) 
(9) 



dV.S. 



dV.S. 



^e^asS {xj* — 1) z(l — z) (z^ + (1 — z)^) 



4e2p^,P_Lj- 



x(167r3 



d^vd^i 



-iq±_ ■ {v — v') r 



2(Tr 



{Pl+e^r 



f J 



(10) 



aemelasS {x^* - 1) 4z^(l - z) 



, ie}P^,P^, 
(Pi + 



x(167r3) 



(27r)e 



'<?-L-(«-i.')2 (^Tr [p*~(w)Z^[+l1'PJ-(t;')Z^W] ) 



Here we have used the identity 

- (Tr {d,U[v)\ V\v') \d,U{v')\ U\v)),^ = gl 



dw+dt;'+ ( Tr 



P'-(v)Z^WtpJ-(i;')Z^[+l 



(11) 



(12) 



where the gauge link iY^+l connects the two coordinate points by means of longitudinal gauge links going to +oo and 
a transverse link at infinity which does not contribute when the appropriate boundary conditions are taken. 

If one integrates over the orientation of Pi, one can replace P±iP±j by I^PjIl This replacement allows us to 
reduce the above expressions into Eqs. (30) and (31) in Ref. '8] which only involve the conventional Weizsacker- 
Williams gluon distribution. 



^ One could also define v = ^{xi + X2) in this process since tlie virtual photon does not have initial interactions with the nucleus target, 
then the respective conjugate momentum is P± = ^{kij_ — k2±) ~ P±. P± is the relative momentum of outgoing partons respect to 
the center of mass frame. Nevertheless, the following derivation remains the same in this case. 

^ In the derivation of Ref. Q , we have employed this as an underlying assumption. 
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Now we are ready to show that the hnearly polarized Weizsacker- Williams gluon distribution can also arise naturally 
in the color dipole model. Since the indices i, j are symmetric, we can decompose the operator expression appearing 
in Eqs ([T0| and (|TT|) into two parts with one part involving only 5ij and the other part being traceless, 



d^vcfv' 
(27r)3 




(13) 



Here a;G'^^-' {x, q±^) is the conventional Weizsacker- Williams gluon distribution while the coefficient of the traceless ten- 
sor xh^^\x,q±) is the so-called linearly polarized partner of the conventional Weizsacker- Williams gluon distribution. 

The physical meaning or interpretation of these two gluon distributions can be better represented in a frame in 
which the two components of the transverse momentum with j — 1,2 or j — x,y are the same. With q^ — q^ in 
this frame, the two symmetric projection operators in Eq. ([T3|) can be written as. 



1i 



1 
1 



1 

2 

1 fO 1 

2 U 



2 \ x'-y 



(14) 
(15) 



where = (IjO) and = (Oj 1) are 2-dimensional unit vectors along a;-axis and y-axis, respectively, which could 
be interpreted as two orthogonal linear polarization vectors for transversely polarized gluons. As shown in Eqs. (jl4p 
and (|15[) . these two symmetric projection operators can also be expressed in terms of the two orthogonal circular 



polarization vectors for transversely polarized gluons, = [Te-i — i e^]/\/2- For the comparison, we also list here the 
antisymmetric projection operator for the polarized gluon helicity distribution. 



1 f Q 



1 

= - i\ 
2 



^x^y 



£ I el — e_e_ 



(16) 



From Eqs. and p^ . it is natural to interpret G*-^^ as a probability distribution to find unpolarized gluons, 
while the polarized gluon helicity distribution could be interpreted as a difference of two probability distributions to 
find positive helicity gluons and negative helicity gluons, respectively. From Eq. dTSl) . it appears that h^^ does not 
have a probability interpretation in terms of the base polarization vectors e^i which are the eigenstates of angular 
momentum operators. □ Instead, it could be interpreted as a transverse spin correlation function to find the gluons 
in the amplitude and complex conjugate amplitude to be in two orthogonal polarization states. In a general frame, 
q^^ = {q^,q\) = (cos sine/)), the projection operator for /i^-* can be written as, 



cos{2(f>) sin(2(^) \ 
sm{2(t)) - cos(20) ) ' 



(17) 



which includes the special case in Eq. (jl5l) when (j) = n/A. Since the projection operator in Eq. (|17p is proportional to a 

rotation matrix of the azimuthal angle, the h^^"^ could also be interpreted as "azimuthal correlated" gluon distributions 
[l3,[I3- Because the gluons in the amplitude and complex conjugate amplitude are in different transverse spin states, 
this kind of gluon distributions could contribute to the observables with cos(2A<?!))-type azimuthal dependence, or 
azimuthal symmetric observables if they come in pairs. 



However, if one chooses different base polarization vectors as e\ = "^(t)t) ^2 ~ ^2^^'"^^' eigenstates of 

angular momentum operators, one can find that Eq. I I15I I becomes ^ (^'^i^i ~ ^2^2^ which would allow us to interpret /i^' as the linearly 
polarized gluon density along the direction of the linear polarization. In a general frame, the polarization vectors are found to be 
= (cos (/>, sin </>) and = (sin <f>, — cos </>) which convert Eq. I|17|l into ^ (^^i ~ ^2^2) ^ well. This indicates that the interpretation of 
the linearly polarized gluon distributions depends on the choice of the polarization vectors. 
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Substitute Eq. dTS]) into Eqs dTO]) and ([TTl), we obtain 



dV.S. 



dV.S. 



aemelas6 (x^,. - 1) z(l - z) (z^ + (1 - zf) [ ^ 

xG^i) (x, gj.) - cos (2A0) x/i^'^ (a;, g^^) 

„ „ „ 8e^Pj 

aeme^as^ (a;^. - 1) 2 (1 - ^) 



(18) 



xG'-^\x,q^) +CO8 (2 A0) x/i^^^ (x, ) 



(19) 



where Ai/) ^ 4>p^^ 4'q± with and 0gj^ being the azimuthal angle of Pi and q±, respectively. This result is in 
complete agreement with the one obtained in Ref. [l^ by using the TMD approach. The coefficient of the cos (2A0) 
term in the above cross section can provide us the direct information of the linearly polarized Weizsacker- Williams 
gluon distribution xh^^\x,q±). It is also easy to see that the xh'^^^ {x,q±) term vanishes if one averages the cross 
section over the orientation of either P± or q± due to the factor cos(2A0). This is transparent when one uses the 
variables P± and q± since they can be interpreted as the relative transverse momentum with respect to the center of 
mass frame of these two outgoing partons and the total transverse momentum of the CM frame, respectively. 

Last but not least, one can see that the contribution from the linearly polarized gluon distribution vanishes if Q = 0, 
i.e., the real photon nucleus scattering only involves the conventional Weizsacker- Williams gluon distribution. This 
is because the real photon cannot generate a cos (2A(^)-type transverse spin correlation that matches the transverse 
spin correlation generated by h^j^\ 

Let us now study the behavior of (a;, q±) in the McLerran-Venugopalan (MV) model[l^ for a large nucleus with 

A nucleons inside. Using the quadrupole results calculated in Ref. [1], one can cast the analytical form of xh^j^\x, q±) 
intoEl 



xh^l\x,q^) = — [S^^- 



2qlq^ 




d^vd'^v' 



dr±r± 



(27r)2(27r)2 



-iq± ■ {v — v' ) 



(Tr [d,U{v)] UHv') [d,U{v')] U\v)),^ 



In- 



' A2 



1 — exp 



(20) 



where J2{qi^r±_) is the Bessel function of the first kind and — asg^Nc^J? In with /i^ — 2§^- For ^ Q 



^^2^ — •- which is identical to xG*-^-' (x, q^) and agrees with the perturbative QCD results. 



we find that xh^^ [x, q±) 

It is important to notice that it scales like A since each nucleon contributes additively in the dilute regime. In this 
regime, the dominant contribution to the gluon distribution comes from a single two-gluon exchange with a transverse 
momentum transfer q^ in the color dipole picture. For the case ^ q\ ^ Q^g one absorbs the In -prj^ factor into 

the definition of the saturation momentum and finds xh''j^\x, q±) ~ ""^qI^" which is an approximate constant. It 

scales like A^^^ since ^ A^^^ as a result of strong nuclear shadowing. It is interesting to note that, in the low 
region, the effect of multiple scatterings between probes and target nuclei can be viewed as or attributed to a 
single scattering with the momentum transfer of order . As compared to the small behavior of the conventional 



N^--l 1 Qi 
m -o 



we find that 



1I 



> 1 



Weizsacker- WiUiams gluon distribution xG'^^^x,q±) 
where we have replaced r± by These gluon distributions obtained in the MV model can be viewed as an initial 



xh^^^ {x,q±) 



In^ln^ 



condition for the small-x evolution. In addition, we can also find that xG^^'^{x,q±) > xh\j {x^q±_) for any value of q±_ 
which ensures the positivity of the total cross section. 



B. Drell-Yan dijet process 

Following the prompt photon-jet correlation calculation in Ref. jS] , it is straightforward to calculate the cross section 
of dijet (97*) production in Drell-Yan processes in pA collisions. The calculation is essentially the same, except for 
the slightly different splitting function since the final state virtual photon, which eventually decays into a di-lepton 
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pair, has a finite invariant mass M. By taking into account the photon invariant mass, the splitting wave functions 
of a quark with longitudinal momentum p+ splitting into a quark and virtual photon pair in transverse coordinate 
space become 



ieMKi{€M\r\)^:;^{Sa-Si3- + (1 - z)Sa+Sj3+), A 

ieMKi{eM\r\)^:;^{5a+Si3+ + (1 - z)Sa-Si3-), A 



(21) 



^tM' ^t,r) = T+ (1 - z)MKo{eM\r\)5^p. 



(22) 



where e\,j = (1 — z)M^, A is the photon polarization, a, (3 are helicities for the incoming and outgoing quarks, and 
z = ^ is the momentum fraction of the incoming quark carried by the photon. 

At the end of the day, for the correlation between the final state virtual photon and quark in pA collisions, we have 



""dp 



2 

t^)^^il-z)z'S^F,^iq^) 




27r2 



PI 




{Pl_ + zq^)2 + 



(23) 



with F,^{q^) = J |l^e-'«^"'^^(TrC/(0)C/t(ri))^^,g^ = fcii + A:2±andPx = (l-2)A:i±-zfc2±. In the MV model, 

with ~ jf'Q% being the quark saturation momentum. qf(xp,fj,) is the integrated 



P-Al^) - ^fee^P (-|t 

quark distribution with flavor / in the proton projectile. Here we used the hybrid factorization which allows us to 
use integrated parton distributions since the proton projectile is considered to be dilute as compared to the nucleus 
target. The flrst term in the curly brackets arises solely from the transverse splitting function in Eq. (PT|) while the 
second term is the sum of contributions from both the transverse and longitudinal splitting functions. We would like 
to emphasize that the above cross section in Eq. (j23p is an exact result regardless of the relative size between q± and 
Pj_. By taking the correlation limit, namely qi_ <^ P±, we arrive at the result which is identical to the one obtained 
from TMD factorization 16] 



""DP 



dyidy2(Pki^<Pk2i 



q±<t:P± 



E 



Xpqf{xp,n)xG^'^\xg,q±) [Hq 



cos (2A0) iJ, 



(24) 



with A0 = - 0,^, xG^^\x,qi_) = §^S^F,^{q^) and 



H. 



asae.m.ej (1 - z) z^ 



PI 



e}{l~z)z^ 2z2 



PI 



2^2 2 p2 



(25) 
(26) 



In this case, the relevant gluon distribution is the so-called dipole gluon distribution as demonstrated in Ref. [tI. [sI. [20| . 



^ To compare with Ref. hJI, one can compute the Mandelstam variables and find that s = {ki + k2)'^ = + — — ^^^"^^ _|_ 
2fcli ■ k2± = u = and t = ^. 
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As discussed in Ref. [l^ , according to the operator definition of dipole type gluon distributions 0, H, HO] , 



(27r)3p 



(27) 
(28) 
(29) 



wfiere tfie gauge link 1 is composed by longitudinal gauge links going to —oo. This shows that the linearly polarized 
partner of the dipole gluon distribution is exactly the same as the dipole gluon distributior0. From Eq. (j29p . with 

the proper normalization, we can also find that the linearly polarized gluon distribution xh^f^ {x,q±) — xG'''^\x, q±). 

Furthermore, one can see that for the prompt photon-jet correlation, the linearly polarized gluon distribution does 
not contribute since vanishes when M = 0. This is also due to the fact that the real photon in the final state 

cannot generate the transverse spin correlation that matches the transverse spin correlation of the incoming gluon in 
the qg — >■ 97 subprocess. It takes two matched transverse spin correlations to get a nonvanish observable effect. 

C. Resummation 

For the purpose of the Collins-Soper-Sterman resummation [2H discussed in Ref. it is also useful to define the 
coordinate expression of the linearly polarized Weizsacker- Williams gluon distribution as follows 



xh1^''ix,b^) 



1 



d qj_e 



-iq±-b± 



2qW^ 



S'^ I xh^l>{x,q^), 



(30) 



and it is straightforward to find that in the MV model 



xh^l^'\x,b^) 



1 



2&Yfei 



62 In ■ 1 



1 — exp 



(31) 



This can be compared to the normal Weizsacker- Williams gluon distribution in b± space defined as zG'^' (a;, 6j 
/ d^qj_e-"}^-''^xG'^^\x,q±), 



x&^\x,b±_) = 



exp 



(32) 



At small b±, xh!"^'''' {x,bj^) is proportional to {^^^ — 26^6-^/6^^ times a constant, whereas xG^^^\x^bi_) behaves as 

In jrgr due to the logarithmic term in Q^g . These properties are consistent with their perturbative behaviors at large 

transverse momentum . 

Similarly for the dipole gluon counterparts, one gets 



In 



In 



WW' 



and 



X&^\x,b,-^- ^'^^ ,2- 



■eM-^Qsqbi]Q 



In 




In 



(33) 



(34) 



Again, in the small b± limit, they behave the same as those Weizsacker- Williams gluon distributions, respectively. It 
is interesting to notice that their large b± behaviors are different. For the dipole gluon distributions, they decrease 
exponentially whereas the Weizsacker- Williams ones have power behaviors. These expressions can be viewed as the 
initial conditions of the resummation discussed in Ref. Ilill . 



There is a factor of 2 between these two distributions in Ref. [l6| due to different normalization. 
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III. SMALL-a^ EVOLUTION OF THE LINEARLY POLARIZED GLUON DISTRIBUTIONS 

In this section, we discuss the smah-a; evohition of the hnearly polarized gluon distributions. We separate the 
discussions into two parts: the first part is on the evolution of the linearly polarized dipole gluon distribution since 
it is trivial and it only involves the dipole amplitude; then we derive the evolution equation for the linearly polarized 
Weizsacker- Williams gluon distribution from the small- a; evolution equation of quadrupoles. 



A. The evolution of the linearly polarized dipole gluon distribution 

According to the definition of the linearly polarized dipole gluon distribution, and the above calculation of the cross 
section of dijct (97*) production in Drell-Yan processes in pA collisions, we know that the linearly polarized partner of 
the dipole gluon distribution is identical to the normal dipole gluon distribution, i.e., xh^^\x,q±) = xG^'^\x,q±^. In 
general, one can write these distributions in terms of the dipole amplitude, namely, the two point function of Wilson 
lines ^ (Tr {U {x {y as follows 

x/if (x,<j^) = xG(^)(a;,g^) = |I#^ / <fx^ f ^e-^-^-"^-)-! (TrC/(x^)C/t(y^))^ . (35) 

27r^Q;s J J {2itY TVc ^ 

The small-x evolution of the dipole amplitude follows the well-known Balitsky-Kovchegov equation [13, [3 which 
reads 



X |(Tr [U{x)U\y)])^ (Tr [C/(a;)t/t(z)] Tr [C/(z)C/t(y)] )^| . (36) 

In the dilute regime, the Balitsky-Kovchegov equation reduces to the famous BFKL equation which leads to the 
exponential growth in terms of the rapidity Y ~ In ^ . 



B. The evolution of the linearly polarized Weizsacker- Williams gluon distribution 



The operator definition of the Weizsacker- Williams gluon distribution can be obtained from the quadrupole corre- 
lator whose initial condition can be throughly calculated in the MV model. In Refs. p2l - [25i |. the small-x evolution 
equation of the quadrupole has been derived and studied analytically. Similarly to the Balitsky-Kochegov equation 
for dipoles, quadrupoles follow BFKL evolution in the dilute regime and reach the saturation regime as a stable 
fixed point. In addition, one expects that quadrupoles should also exhibit the same geometrical scaling behavior as 
dipoles. Recently, using the JIMWLK renormalization equation [1^, H^l, the first numerical studies of the small- a; 
evolution of quadrupoles indeed observe evidence of traveling wave solutions and geometric scaling for the quadrupole. 
According to and Ref . 0, H, US] : the Weizsacker- Williams gluon distribution can be written as 

xG%^{x,k^) = -^J (02 ^ e-'=-(''-"') (Tr [d^U{v)\ U^v') [dW{v')\ U\v))^ . (37) 
= ^xG(i)(a;,g^) + \ (^2^ S^^^ xh^l\x,q^). (38) 

The evolution equation for the correlator (Tr [9W(u)] ^{v') [d^U{v')\ W{v)')y can be obtained from the evolution 
equation of the quadrupole (Tr (C/(a;i)C/^(a:']^)/7(a;2)C/^(a;2)))-j^ by differentiating with respect to x\ and Xj, and 
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then setting x\ — X2 = and = x'l 



Then the resulting evolution equation become^ 



d 

— ^ir \a^v \v)\ v\v^) V [v^ )\V^\v),^ 



(Tr [d'U{v)] U^{v') [d^U(v')] U^v))^ 

— {TT[d^Uiv)] uHv') [dWiv')] uHv)), 
V — v'Y {v — zY 



27r2 



(fz. 



d z±_ 



{v — zY{z — v'Y 
1 {v-v'f 



27r2 J ^ Nc{v- zf{z-v'Y 
X {(Tr [U(v)U\v') [d'U{v')] U\z) 



asNc 
' 27r2 



Sz±_ 



{v - v') 



l\2 



{v~ zY{z-v'f 



V — v'Y [v — zY'' 
Tr {U{z)U\vY\)y - (Tr \U{z)U\v') {&U{v'Y\ U\vY\ Tr \U (v)U\zY\) y} 
v' - vY {v' - zY 



u)2 [v' — zY 



X {(Tr [[d'Uiv)] U\z)U{v')U\v)] Tr [Uiz)U\v')])y - (Tr [[d'U{v)] U^v')U{z)U^v)] Tr [U {v' )U^z)]) y} 



' 47r2 



{v-v'Y 



zY{z-v'Y 



X {(Tr \U{v')V\zY\ Tr \U{z)V\v'Y\ )y + (Tr \U{v)U\zY\ Tr \^{z)V\vY\ )y 

- (Tr {U{v')U\vY\ Tr {U(v)U\v'Y\)y - . (39) 
The evolution equation of the Weizsacker- Williams gluon distributions can be obtained by contracting the above 



correlator with Sij and the one for the linearly polarized partner by contracting with ( 2 



^ . Although the 



expression is quite complicated in general, the result gets simplified in the dilute regime as in Ref. [23]. In the dilute 
regime, the correlator which yields the Weizsacker- Williams gluon distribution can be reduced to a simple form in 
terms of r(w, w') 



(Tr [d'U{v)] uHv') [d^U{v')] UHv))y = ^9ia^,r(z;, z;')y , 
where ^f^T{v, v') is the leading order dipole amplitude which satisfies the BFKL equation 



— r(a;i,a;2)Y 



d^z- 



{xi 



X2) 



[r{xi, z)y + r(z, X2)y - T{xi,X2)y] 



(40) 



(41) 



27r2 J (xi - Z)2(X2 - Z)2 

In the dilute regime where the gluon density is low, we know that the Weizsacker- Williams gluon distributions, 

xG^^-* (a;, (jx) and xh'\i'' {x,q±), as well as the dipole gluon distributions all reduce to the same leading twist result. 
Therefore, despite of the distinct behavior in the saturation regime, we find that all these four types of gluon distribu- 
tions follow the BFKL equation in the dilute regime where the gluon density is low. The physical consequence of this 
results is that the linearly polarized gluon distributions should be as important as the normal gluon distributions in 
the low- a; region since they also receive the exponential rise in rapidity Y due to the BFKL evolution. Furthermore, 
according to the discussion in Ref. [29, 30], the BFKL evolution together with a saturation boundary can give rise 
to the geometrical scaling behavior [3ll - l33| of the dipole gluon distribution. Since the quadrupole evolution equation 
also contains the same property as discussed in Ref. [23, [2^, the Weizsacker- Williams gluon distribution and its 
linearly polarized partner should exhibit geometrical scaling behavior as well, although their evolution equations are 
much more complicated in the saturation regime. In terms of the traveling wave picture (29l . [s^ j for the evolution 
of dipoles and quadrupoles, the velocities of the traveling waves for dipoles and quadruples are identical, since the 
velocity is determined by BFKL evolution. This implies that the energy dependence of the saturation momentum 
~ Qq[xo/x)^ with Qo = IGeV, xq = 3 x 10^^ and A = 0.29, should be universal for all these four different gluon 
distributions. 



IV. CONCLUSION 



We perform the color dipole model calculation of the cross section of DIS dijet and Drell-Yan dijet processes, and 
demonstrate that the linearly polarized partners of the Weizsacker- Williams and dipole gluon distributions naturally 



This evolution equation involves derivatives of the Wilson lines and complicated kernels which make it very hard to solve directly. 
However, one can extract the evolution information by numerically solving the evolution equation for quadrupoles first and then making 
numerical differentiation and identification of coordinates. 
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arise in these processes. This result is in complete agreement with Ref. jl^ and implies that the measurement of the 
cos (2A0) asymmetries in these dijet processes can be a direct probe of these two different linearly polarized gluon 
distributions. In addition, the small- a; evolution studies of the linearly polarized gluon distributions reveals that they 
also rise exponentially as function of the rapidity at high energy and they should also exhibit the geometrical scaling 
behavior as the normal unpolarized gluon distributions do. 
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